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Abstract 

The properties of baryon resonances are extracted from a complicated process of fitting sophis- 

CN ' ticated, empirical models to data. The reliability of this process comes from the quality of data 

> \ 

■ and the robustness of the models employed. With the large of amount of data coming from re- 



cent experiments, this is an excellent time for a study of the model dependence of this extraction 
process. A test case is chosen where many theoretical details of the model are required, the Sn 
partial wave. The properties of the two lowest A^* resonances in this partial wave are determined 
using various models of the resonant and non-resonant amplitudes. 
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I. INTRODUCTION 



The problem of extracting baryon resonance parameters from observables measured in 
scattering experiments is of fundamental importance. In order to understand the relevant 
degrees of freedom at low energy and their interactions in the few dozen A^* states found to 
date, reliable and objective information on their properties extracted from data is required. 
Models of the baryon spectrum are usually compared to masses defined in terms of effective 
Breit-Wigner parameters. A more sensitive test of such models is the comparison of total 
widths and branching fractions to open channels of each resonance to those extracted from 
data. These quantities are also defined in terms of effective Breit-Wigner parameters, but 
in general show more sensitivity to the method used to fit the data. What is not commonly 
appreciated is the dependence of these resonance parameters on the method used to extract 
them from a fit to the data. For example, it is not uncommon to see such parameters referred 
to as "data" in the literature. 

The primary goal of this paper is to examine the model dependence in this process in 
a carefully chosen test case, states in the Su partial wave in vrA^. States in this partial 
wave have both significant physics interest and significant uncertainty in their parameters 
as reported by the Particle Data Group (PDG) in their Review of Particle Properties jlol |. 
Most of the literature considers the lowest energy resonance in this partial wave [5*11(1535)] 
as a three-quark state within the quark model ^Ij]. However, the substantial branching 
fraction of 5*11 (1535) to r]N has, given the small phase space available for this decay, led 
to the interpretation of this state as a meson-baryon molecule P]. Lattice calculations of 



the masses o: 
Ref. H and 



these states have recently been developed ^, J|. The calculations of both 



are quenched, and find masses for 5*11(1535) similar to the values fit from 
experiment [lO']. In addition, constituent quark models which describe these 5*ii resonances 
as (P-wave) orbital excitations with = l/2~ differ in their prediction of the amount 

n 

of mixing between the two possible quark-spin states pj, because of different models of 
the short distance interactions between the quarks. A comparison between partial widths 
calculated in a given model and those extracted from data can be used to determine the 
nature of these short distance interactions. This makes the dependence of these partial 
widths on the method used to extract them from fits to data important, and relevant. 
These important microscopic issues are unable to be settled with the uncertainty in the 
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full and partial widths presently reported by the PDG, which vary widely among various 
studies. At one extreme, analyses of the ttN elastic scattering data tend to give a full width 
of about 120 MeV for 5*11(1535), with vrA^ as the dominant decay branch ^,Q|. At the other 
extreme, threshold eta photo-production data gives a much larger full width F ~ 210 MeV, 
with TjN as the larger decay branch jSi]. Only a unified treatment of all the data can provide 
a consistent picture. 

The purpose of this paper is to examine the model dependence of 6*11 resonance param- 
eters by extracting them from several fits to a single set of partial wave amplitudes for 

0n 
and vrA^ r]N [9j scattering in the 5*11 partial wave. These fits have differ- 
ent levels of sophistication in their of description of the scattering matrix. The analysis of 
these channels in this partial wave is chosen because the proximity of two overlapping reso- 
nances [5*11(1535) and 5ii(1650)] and the rjN channel threshold makes it an interesting and 
non-trivial example. Differing treatments of the re-scattering processes are a key difference 
among the models studied here and these effects should be most important near a threshold. 
Although the number of open channels is much larger than two, the two channels that have 
been chosen are the most important and account for 90% (1535) and about 80% (1650) of 
the decay width when all channels are included [lol |. 

Breit-Wigner models can get good fits to single channel data (e.g. vrA^ rjN jll|. 
7A^ — > T^N [12]), but have non-trivial uncertainties and at times surprising results. In 
simplifying an intrinsically multichannel problem into a single channel model, information 
from other reactions must be included with accompanying uncertainties. In addition, there is 
no single commonly accepted form for the resonance shape. When a somewhat large (~10) 
set of asymptotic states can couple to each A^* with high probability, including a large 
sample of reactions and application of the unitarity constraint are extremely important. 
Common multichannel techniques include if-matrix 14 1 and Carnegie-Mellon Berkeley 
(CMB) models. Each has an established way of including a variety of interactions 

respecting multichannel unitarity with varying ability to include dynamics. 

Three methods are used here to perform the extraction of resonance parameters via a 
fit to the par tial wave amplitudes. The most constrained is the CMB model, which was 
developed [16[ to describe i^N elastic and inelastic scattering, and extended and modernized 
in Ref. j^. This approach has multi-channel unitarity, and the scattering matrix (a matrix 
in the space of channels) has the required analytic structure. Resonances are modeled as 
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"bare" poles, which are "dressed" by couphng to the open asymptotic channels. 

The second, simpler approach is to describe the unitary T matrix using a real symmetric 
K matrix, which is in turn written directly in terms of resonance parameters. This approach 
maintains multi-channel unitarity, but does not satisfy analyticity constraints. The third, 
simplest approach is to build the T matrix directly from a sum of relativistic Breit-Wigner 
forms for the resonances. In this case, minimal constraints on the scattering amplitudes are 
available. Below tttt production threshold, the overall amplitude can be made unitary using 
Watson's theorem. However, the rjN channel must be included for a good description of 
the 5*11 partial wave, 

Breit-Wigner models used here are neither unitary nor analytic. 

Non-resonant amplitudes must be added to account for scattering processes which do not 
involve s-channel N* resonances. Although a fairly small set of diagrams describing non- 
resonant processes can be identified at lower total energy {W ~ 1.3 GeV), the set of possible 
diagrams grows rapidly as the total energy increases to 2.0 GeV. To date, no published 
work has included all the relevant diagrams. Empirical descriptions of the non-resonant 
amplitudes can be chosen since each resonance has a strong signal in at least one of the 
reactions studied, and most publications find the non-resonant amplitudes to be smaller 
than the resonant amplitudes. In addition, the most basic physics assumption is that the 
resonances come from the long-distance part of the interaction and are seen in the sharp 
energy features of the data, while the non-resonant amplitudes arise from the short-range 
part of the interaction and provide smooth energy dependence. This implies that, in this 
case, the influence of the choice of non-resonant amplitudes on the extracted resonance 
properties should be small. The "distant poles" model of the non-resonant amplitudes is 
designed for use with the CMB model. Bare poles well below and well above the thresholds 
for the channels being studied, i.e. distant poles, are fit to the partial-wave data with 
methods very similar to those used for the resonance poles. In the CMB model they provide 
a smooth, analytic, unitary T matrix in the physical region; they can also be incorporated 
into a i^-matrix model to provide a unitary scattering matrix. While the poles and cuts 
in the left-hand part of the complex s-plane required by crossing symmetry are difficult to 
fully implement in the CMB model, their effects in the physical region can be simulated by 
a particular choice of strong form factor, and by the use of sub-threshold poles. Another 
method adds to the resonant amplitudes polynomials quadratic in the parameters Xc = 
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{s — Sth.c)^*^^ for each exit channel c, which can yield a unitary T matrix for the i^'- matrix 
approach. 

The procedure used here is to calculate, using the CMB, i^-matrix, and Breit-Wigner 
approaches, the matrix Tab for scattering between the channels {a, b} G {ti'N, rjN} in terms 
of a set of parameters describing the resonant and non-resonant amplitudes. Different treat- 
ments of the background are explored, giving a total of eight different models. A statistic 
is minimized by varying theseparameters to best fit the partial wave amplitudes for the 
vrA^ — s> TiN 0] and vrA^ — > rjN js] reactions. 

The goals of this work are somewhat limited in order to make a clear statement about 
model dependence, and for this reason the resonance parameters which result should not be 
considered for use in understanding the microscopic structure of the resonance states. 
The channel space is greatly simplified to only the two principal channels, leaving out other 
channels such as ttA that are a much smaller part of the amplitudes, but which must be 
included for the optimal values. From ten-channel fits to both ttN and reactions, the 
5*11(1650) width is known to be ~50% larger than in a two-channel fit. In addition, the 
sophisticated i^-matrix techniques developed by the Giessen group [l3| are not used here. 

II. MODELS OF THE SCATTERING MATRIX 

In this section the form of the scattering matrix in the various models is described, 
along with the implementation of the two forms of background described above. All of the 
models used here require partial wave amplitudes (PWA) as input. These extract the energy 
dependence of amplitudes with specific isospin, parity and angular momentum (e.g. Sn) 
from the large number of experimental data points. About 20 years ago, significant efforts 
were made to find model independent methods for partial wave analysis of irN elastic data [3| ■ 
More recent efforts by the George Washington University (GWU) group incorporate most 
of the theoretical constraints developed previously. Although the fits in this work do not 
use the most recent GWU work, results using the most recent PWA would not produce 
different conclusions. Since the number of data points for the tcN rjN reaction is much 
smaller, there is more model dependence. These amplitudes are discussed in Ref. jol, using a 
method that that accurately couples the nN elastic and inelastic data with minimal model 
dependence. Until another study is done, this is the only PWA result available for vrA^ rjN. 
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A. CMB model 



In this work the 5*11 partial wave in vrA^ elastic and iiN — > rjN scattering is described 
using only the two light resonances 5*11(1535) and 5*11(1650). In the CMB model background 
{t- and M-channel) processes are simulated using one high-energy and two sub-threshold s- 
channel poles, which are treated identically to the resonances in order to preserve the analytic 
structure of the T matrix. The transition amplitudes 0, Q| are 

5 5 

= 5Z 5Z fa{s)\/ Pa{s)-iaiGij{s)-ibk\/ Ph{s)fb{s) (l) 
i=l j=l 

where {a, G {ttA^, r^A^} label the channels, and i,j label the poles, which can represent 
the resonances or the background. The model used in this paper has five poles (two for 
real resonances and three for background resonances). It is easy to add channels with this 
formalism if the input data is available. The phase space factor Pa{s) has the form 

Pais) = ^, (2) 

where pa is the cm. frame momentum and jai is the real-valued coupling constant of 
resonance i to the channel a. The form factor fa in Eq. (Q) is unity in the 5*-wave. 

The dressed propagator Gij{s) allows resonance i to couple to resonance j through re- 
scattering, and is the solution of the Dyson equation 

G,,{s) = G°.(s) + Gl{s)J:,i{s)G,,{s). (3) 

The bare propagator G° (s) has the form 

GUs) = (4) 
soi- s 

where = +1 for the resonances and the high-energy background pole, and soi is the bare 
mass of the resonance (background pole). One sub-threshold pole has Cj = -|-1 (repulsive), 
and the other has Cj = — 1 (attractive). 

The self energy T,ki in Eq. ^ is the sum over channels c 

2 

Sfc/ = ^7cfc0c(s)7d, (5) 

c=l 

where 0c(s) is the channel propagator, which plays a central role in the CMB model. The 
sum is over the two channels in this simplified problem. Two-pion channels such as pN 



can be included by treating them as quasi-two-body channels, which results in an enlarged 
T matrix. The re-scattering sum in Eq. will then extend over the additional channels. 
The imaginary part of 0c is the product of phase space, Pc{s), with the square of the form 
factor for channel c (this quantity will be seen in the other models below) and the real part 
is obtained through a once-subtracted dispersion integral jof Q • The desired analyticity of 
the entire amplitude is then assured. 

Unitarity comes about because of the democracy of the re-scattering in the Dyson equa- 
tion, Eq. (jni). In the multichannel context, unitarity is expressed through a generalized 
optical theorem, 

lmT,, = J2T:jci> (6) 

c 

Unitarity is also related to the properties in the complex plane through the discontinuity in 
the amplitude across the right-hand cut. 

The dressed propagator Gij{s) can be found by solving algebraically the matrix equation 
Eq. (ini), with the result 

H,,{s) ^ [G-\s)],, = [G\s)-% - = '-^^6., - (7) 

and the matrix G{s) can be found by inverting H{s). This completes the necessary ingredi- 
ents for calculating the T matrix in the CMB model of Eq. in terms of three parameters 
for each resonance or background pole. These are the bare mass squared so,i, and the cou- 
phng strengths •y^N,! and 7r,Ar,i. Since there are five poles, this model has a total of fifteen 
parameters. 

Once this T matrix is fitted to the partial wave data by varying these fifteen parameters, 
baryon resonance parameters are extracted from the resonant part of the T matrix. The 
procedure for doing this is described in what follows. 



B. Extraction of resonance parameters in the CMB model 

Poles in the resonant part of the T matrix occur at complex values of s where the denom- 
inator of the resonance propagator vanishes. A search program finds these pole locations, 
and these are the model-free output (see Ref. 0] for details). This is possible because the 
amplitude has reasonable properties for complex values of s. 
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However, it is important to derive effective Breit-Wigner parameters for each resonance 
from this modeL Once the pole positions are found, the matrix if(spoio) is diagonahzed 
to ehminate interference between resonances. The denominator of the resonant part of the 
T matrix in the vicinity of each pole is then expanded in s. The constant term becomes 
the dispersive correction to the mass, and the term linear in s becomes the width of a 
generalized relativistic Breit-Wigner form. Thus, although the energy dependence of a 
resonance in the CMB model is much more sophisticated than the normal Breit-Wigner 
shape, all characteristics of a resonance can be expressed in this commonly used form. 



C. CMB model with polynomial background 

This is a modified CMB model which uses a polynomial function to describe the back- 
ground. Care must be taken to maintain unitarity of the T matrix. Using a technique from 
Ref. jl7|, the full S matrix is written 

S = B^RB, (8) 

where R and B are the resonant and non-resonant S matrices, respectively. Here R = I+2iT 
uses the CMB-model T matrix, and i? is a unitary matrix that need not be symmetric. The 
background matrix B is constructed from a real symmetric matrix 

E= + (9) 

where for the two channel case of interest 

K = 

rp I ^2 1^ ^ I J^/ ^2 

or 

^ab — ^abXb + ^afe^65 (H) 

Xb = \J s — {rriN + vrthf' for s above exit channel threshold and zero otherwise, and the 
coefficients kab and A;^^ are real. Note that the off-diagonal terms in are zero for s 
below TjN threshold. This ensures that the diagonal and off-diagonal elements of R are 
not mixed, so the vrA^ rjN amplitude vanishes below riN threshold. This model has six 
real parameters describing the background, and six describing the resonances, for a total of 
twelve real parameters. 
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(10) 



Once the 5* matrix in Eq. (jS} is formed, it can be easily converted using S = I + 2iT to a 
T matrix that can be fit to the partial wave amplitudes by varying these twelve parameters 
in the usual way, and the extraction of the baryon parameters is identical to that of the 
original CMB model. 



D. K-Matrix models 



A unitary 5* matrix can be constructed from a real, symmetric X-matrix via T = K{1 — 



mm 

Is y 1 



for investigating 



iK)~^. This method is very common in the literature 

hadrons in reactions. Recently, effective Lagrangian models |13[ have been used with a K- 
matrix for studies of N* states. Non-relativistic reductions of Feynman diagrams make these 
studies more sophisticated than what is presented here. The earlier work of Moorhouse and 
collaborators [l9| is very similar to the method used here. 

The differences between the CMB and fC-matrix models are well known. A common 
problem with all ii'-matrix methods is the difficulty in maintaining analyticity. Only the 
work of Longacre [l^ accomplished this. Thus, pole positions and residues can be obtained 
from the CMB fits and not from the i^'-matrix fits. Many researchers feel these values have 
less model dependence. Any microscopic model (e.g. lattice QCD) which has the analytic 
properties as an output can only use the CMB results. While the CMB model has a full 
treatment of re-scattering, only the imaginary part of the channel propagator is included 
in the ii'-matrix model. The choice of the Sn resonances near rjN threshold as a test case 
is motivated by the expectation that these re-scattering effects are most important close 
to threshold. As discussed above, the i^-matrix method can be incorporated directly into 
an effective Lagrangian formulation 2. As a result, the non-resonant amplitudes can be 
added diagram by diagram. This is a significant advantage over the present version of the 
CMB model, which is formulated in terms of arnplitudes rather diagrams. In this work, we 
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use the amplitude form of the i^-matrix model 

Although the CMB model is preferred on theoretical grounds, the i^-matrix model has 
practical advantages. It is more commonly used, and the less complicated treatment of 
interference and re-scattering effects simplifies the fits to partial-wave amplitudes. 

As in the CMB model, two sub-threshold poles and one high-energy pole can be used 
to describe the non-resonant part of the scattering matrix. The contributions from the 
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resonance and non-resonant poles are 



r.i , N V Pa{s)fa{s)laahih{s) ^ Pbjs) , . 

[Si - S) Ci 

Here the form is similar on the surface to what is used in CMB model. The phase space Pa{s) 
is given by Eq. fa{s) is the form factor (none is needed for S-wave), jai is the coupling 
constant of the resonance or background pole i to channel a (but with units different from 
the CMB model), and Si is the position of the i-th (real) pole. The signs Cj are the same as 
those used in CMB model, with e = +1 for the resonances, the high-energy background pole, 
and one sub-threshold background pole, while the second sub-threshold background pole has 
e = — 1. The ii"- matrix is formed from a simple sum over the resonance and background 
poles 

Kabis) = J2 ^abis) (13) 

iGpolcs 

When the i^'-matrix is converted to a T-matrix, re-scattering ( "dressing" of the resonance) 
comes about naturally and the resonance gains a finite width. In this model the ii'-matrix 
is written directly in terms of the resonance parameters, so that Mj = ^/sl is the mass of the 
resonance associated with the pole at Sj, the partial width of the i-th resonance to channel 
a is 

T. Pa{Si)fa 

and the total width is the sum of the partial widths. As discussed above, the differences 
between CMB and i^'-matrix formalism are seen only when interactions with other reso- 
nances and non-resonant amplitudes occur. These differences are expected to be maximized 
close to a channel threshold, which is the reason for the choice of the Su partial wave with 
resonances near rjN threshold. In this case, the CMB model is expected to give more reliable 
results. 

An alternate description of the background is using the polynomial form of Eq. ()lUp. 
which is added to the resonance K matrix 

K,,= Yl Kbis) + Kg. (15) 

i€ resonances 

With either method of treating the background, these ii'-matrix models have the same total 
number of parameters as the CMB models, 15 for the distant-poles background and 12 for 
the polynomial background. 
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E. Breit-Wigner models 



There are two kinds of Breit-Wigner model used here along with two kinds of non-resonant 
amplitude parameterization. In each model, the resonant and non-resonant T-matrices are 
summed, 

Tab{s) = T'^b,resis) +Tab,nonres{s), (16) 

iG resonances 

This form is inconsistent with unitarity without special effort. 
The resonant amplitude T^bresi^) is given by 



T^,res(^) = , e^^f. (17) 

ei [Si - s) - ^^ysTi 

The signs Cj are the same as in the CMB and i^-matrix models, Sj is the pole position (with 
Mi = ^/si for resonance poles), and Fj = J2a is the total width of resonance i. Using a 
procedure similar to that of Ref. jl^, an arbitrary phase exp{i6f^) is associated with each 
pole in each channel. As shown below, the freedom to fit these additional phases is crucial 
to achieving a good fit with a Breit-Wigner model to the Su partial-wave T matrix elements 
in these channels. 

The width commonly used to describe a non-relativistic Breit-Wigner form is the energy- 
dependent width 

Kns) = TUs.)^y (18) 

with the phase-space factor Pa{s) as in Eq. Q. (Note that this is not the same as the 
conventional non-relativistic Breit-Wigner energy dependence.) An alternative, relativistic 
form results from the assumption that the numerator of the contribution of resonance i 
to the i^'-matrix in Eq. (jl2|) is the same as that in the Breit-Wigner form in Eq. (jl7|) . If 
the partial width used in the if-matrix formalism, Eq. (jl4p . is generalized to an energy- 
dependent function by replacing Mj with ^/s, then the i^-matrix numerator of Eq. (|T^ is 

Vpais)fa{s)'yaabiVpb{s)fb{s) = \l T ai^/s sj T u^fs = \^r~Y\iy/s, (19) 

which is the same as the numerator in Eq. (fT7j) . This suggests the use of 

r-'(.) = r..(..)^f^, (20) 

V* Pa\Si) 

for the energy-dependent partial width in a relativistic Breit-Wigner form. In what follows 
models using both r^f (s), labeled BWn.r., and FJjf (s), labeled BWrei, are used in fitting the 
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Breit-Wigner form in Eq. ()17|) to the 5*11 partial-wave T matrix elements, and the results 
are compared. The relativistic form is the same as that advocated by Chung et al. I^i^. 

The first non-resonant form uses distant poles to describe the background contributions 
to the T matrix, so that the total T-matrix is a simple sum 

= n,{s), (21) 

iGpolcs 

This form of background has two low energy poles and one high energy pole, the same as was 
used with the CMB and i^'-matrix models. To get a good fit, each of these terms requires 
an additional phase in each channel, as described above. With the background described 
by distant poles, there are 15 parameters (mass, width, and vrA^ branching fraction for each 
term) associated with the five poles, and ten phases, for a total of 25 real parameters. Of 
the models used in this work, this model has the largest number of parameters. 

As with CMB and i^-matrix models, it is possible to describe the non-resonant back- 
ground by a polynomial function, so that 

Tab{s) = T'^b,res{s) +Tab,nonves{s), (22) 

igresonances 

where the background T-matrix elements are the polynomials 

^a6, nonres ~ f^ab^b ~l~ ^ab^^b' (^^) 



As with the other models, Xb = y s — {niN + rnhf for s above exit channel threshold and 
zero otherwise, but now the coefficients tiab and are complex. In this model there are 
six resonance parameters associated with the two poles, and 8 background parameters (2 
complex numbers for each of two channels), for a total of 18 real parameters. 



III. RESULTS 



In order to understand the complications in the Su partial wave, model predictions for 
a hypothetical isolated resonance are first considered. Figure Q compares the results of all 
four models for the scattering amplitude of an isolated single-channel S'-wave resonance. 
An S'-wave resonance is chosen with a mass of 1710 MeV and width of 215 MeV in the 
vrA^ — s> rjN reaction, which avoids any complications with centrifugal barriers. Agreement 
among the different formulations is striking, despite the deliberate choice of a resonance 
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1.4 1.6 1.8 1.4 1.6 1.8 2 

Energy (GeV) 

FIG. 1: Real parts (left panel) and imaginary parts (right panel) of vrA^ rjN resonant scattering 
amplitudes calculated with an isolated S-wave resonance of mass 1710 MeV and width 215 MeV, 
in each of the four models used here. 

mass close to rjN threshold. The i^-matrix and CMB models have identical forms for an 
isolated resonance. The relativistic Breit-Wigner is chosen to be identical with the CMB 
model in this limit. Even the Breit-Wigner amplitude with the non- relativistic width is 
in good agreement with the other other forms. If the complex amplitudes are plotted on 
an Argand diagram (with the real part of the amplitude on the horizontal axis, and the 
imaginary part on the vertical axis), all will show the typical counter-clockwise circular 
motion of a resonance. Non-resonant background will cause a shift or a distortion of this 
basic shape, and inelasticity will decrease the radius. Interference between resonances can 
have significant effects. 

A partial wave amplitude with this W dependence is the most visible signature of a 
resonance, producing a peak in the total cross section with an appropriate width. Such a 
peak is seen in the total cross section for ttN — > rjN, and is often interpreted as a resonance. 
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1.6 1.8 1 1.2 

Energy (GeV) 

FIG. 2: Real parts (left panels) and imaginary parts (right panels) of ttN ttN (upper panels) 
and ttN r]N (lower panels) scattering amplitudes, T, in the Sn partial wave. The curves in the 
lower panels are the single-channel, one-resonance Breit-Wigner fit to these amplitudes described 
in the text. Dotted lines give the resonant, dashed lines the non-resonant, and solid lines the total 
scattering amplitude. Partial-wave amplitudes are shown with error bars. 



The 5*11 partial wave amplitude dominates the total cross section for this reaction, and is 
shown in the bottom two panels of Fig. El Since the higher mass resonance S'ii(1650) couples 
weakly to this channel, this appears to be an isolated resonance with small non-resonant 
amplitudes. 

Single-channel fits of the vrA^ —>■ rjN amplitudes have been made, which are shown as 
curves in the bottom two panels of Fig. |21 The results are similar to those of single- channel 
fits with a^ Breit-Wigner energy dependence to the total cross sections for Tr^p rjn and 
7P — > rjp [8]. The data can be fit using an S-wave Breit-Wigner form if a small, but important 
contribution from non-resonant amplitudes is included and the energy region used in the fit 
is truncated. This rough fit can be used to determine the product of the initial and final 
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state couplings to the resonance, requiring information from more complete fits to account 
for the missing decay strength of the S'ii(1535) resonance to other channels. The single- 
channel, one resonance fit shown in Fig. |21is made with a Breit-Wigner resonant shape with 
the relativistic form of the width, the polynomial form of the background, and a single phase 
multiplying the resonant amplitude, truncating the fit at 1700 MeV. The result is a mass of 
1543 MeV, a width of 134 MeV, and a vrA^ branching fraction of 48%. This shows that it is 
possible to find fit parameters that roughly agree with the more complete models, but these 
results are of uncertain value given the strong model dependence. 

The elastic tiN amplitude (the upper two panels of Fig. E]) is more complicated, showing 
two structures overlapping in energy. The resonance S'ii(1650) couples strongly to vrA^ with 
peaks at ~ 1.7 GeV in Im(T) and at W ~ 1.65 GeV in Re(T). Although this rapid energy 
dependence is a clear resonant signal, there is also a nontrivial non-resonant amplitude which 
greatly distorts the typical Argand diagram. The energy dependence of 5*11(1535) is more 
complicated in irN elastic scattering because of the strong coupling to the riN channel at 
its threshold (ly =1.487 GeV), which is within the resonant shape. This produces a cusp in 
the real part at threshold in addition to a peak in the imaginary part at approximately the 
resonance mass. Given this complicated structure, analyticity of the scattering amplitude 
and dispersive scattering effects can be expected to be most important for this state. For 
this reason, the elastic amplitude includes many of the important dynamical effects studied 
here. 

After this somewhat pedagogical introduction, the main results of this work are now 
presented. The fits to the partial-wave data for the T-matrix elements T.„N^nN and 
Tt^n^tiN for each resonant-non-resonant model are shown in Figures El to ^| Unlike most 
previous papers studying A^* resonances, which show only the full amplitude, the non- 
resonant amplitude (with all resonant couplings set to zero) and the resonant amplitude 
(with all non-resonant couplings set to zero) are also shown here. The effective Breit-Wigner 
properties of the two resonances extracted from these fits are given in Tables HI and ITTl Errors 
in the first four columns of results are determined from the fitting uncertainties only. 

A common criterion in fits is the value of given in the last row of the tables. A 
comparison of the two-channel fits shows that the lowest values of are attained using the 
Breit-Wigner models with distant poles background, the more theoretically sophisticated 
models are in between, and the Breit-Wigner models with polynomial background have the 
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model 


CMB 


K 




BW,el 


CMB 


CMB 


PDG 






matrix 






all ttN 


all ttN, 




5'ii(1535) 
















Mass (MeV) 


1539±1 


1533±1 


1549±2 


1558±1 


1547±3 


1539±5 


1520-1555(1535) 


Width (MeV) 


135±4 


115±3 


142±9 


143±4 


131±19 


122±20 


100-200(150) 




29±1% 


34±1% 


67±5% 


67±1% 


34±4% 


39±5% 


35-55% 


S'ii(1650) 
















Mass (MeV) 


1682±1 


1685±2 


1648±5 


1637±2 


1690±12 


1684±15 


1640-1680(1650) 


Width (MeV) 


144±3 


190±5 


147±10 


145±4 


227±40 


227±58 


145-190(150) 


BnN 


80±1% 


77 ±1% 


74±5% 


79±1% 


75±3% 


75±3 


55-90% 




3.8 


3.7 


1.5 


1.9 


3.6 


5.6 





TABLE I: Results for resonance parameters from fits to the T matrix elements for vrA^ ttN, r]N 
in the partial wave, using the CMB, i^-matrix, and Breit-Wigner (BW) models described in 
the text. The third last column shows the results of a fit including ttN ttN, 7]N , pN, vrA, aN 



and ttN* partial wave amplitudes |9|]. The next to last column shows results for a fit of the ttN 
data and jN irN, r]N data 12 ill, and the last column shows the range in the central values and 
estimated values from the PDG 10|. Non-resonant contributions to the T matrix are described in 
terms of distant poles in all cases. Since the primary results of this work come from a two-channel 
model, only the branching fraction to ttN is given, since BrjN = 1 — BnN- The last row gives the 
per data point of each fit. 



highest values. Breit-Wigner models commonly have better fits because they are most 
often applied to study single-channel reactions, but are much more likely to have additional 
local minima close to the global minimum, reflecting the lack of theoretical constraints. The 
extra parameters required to get a good fit unfortunately obscure the physics results. The 
extraction of physically meaningful results for resonances depends more on the quality of the 
theoretical constraints placed on the fit than on the quality of the fit itself. In some cases 
the fit function is not as sharp as the data, e.g., the CMB and i^'-matrix model lm{TnN-*nN) 
amplitudes with polynomial background, although this is not apparent from the total 
which is summed over all data points. All models are able to match the shape of the cusp 
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model 


CMB 


K matrix 


BW„.r. 


BWrel 


511(1535) 


Mass (MeV) 


1526±2 


1533±1 


1539±2 


1538±2 




Width (MeV) 


112±6 


119±3 


130±6 


130±6 






30±2% 


33±1% 


39±1% 


38±1% 


5ii(1650) 


Mass (MeV) 


1688±2 


1682±2 


1648±2 


1647±2 




Width (MeV) 


193±6 


184±5 


109±5 


109±5 






78±2% 


75±1% 


51±1% 


51±1% 






5.0 


3.9 


5.0 


5.0 



TABLE 11: Caption as in Tabled except non-resonant contributions to the T matrix are described 
in terms of polynomials. 

at at TjN threshold in Re{TT,N^T,N), since this is a property of interfering amphtudes and 
there are many possible solutions. 

In the figures more detail is shown than is customary. For all models, the amplitudes 
are separated into resonant and non-resonant parts. The resonant amplitudes provide rapid 
energy fluctuations and a rough match to the data, especially for ttN —>■ r]N. The non- 
resonant amplitudes are generally smooth. However, analyticity constraints require a cusp 
at T]N threshold, a feature of all non-resonant amplitudes using distant poles background. At 
first glance, the fits all look similar. More careful inspection reveals differences in the detailed 
balance between resonant and non-resonant amplitudes and in the channel coupling effects. 
One of the most striking features is seen in the vrA^ — > vrA^ amplitude for the relativistic 
Breit-Wigner model with distant poles background. Fig. IHl The imaginary part of the non- 
resonant amplitude is much larger at ~ 1.7 GeV due to coupling to rjN than in the other 
models, a coupled channel effect. The resonant part of the amplitude is then smaller than 
in the other models and this model has a very small for 5*11(1650). 

The resonance parameters extracted using the various models are compared in Tables H] 
and |n] Variations among the results for the resonance masses, the full width and the 
vrA^ branching fraction are all significant given the estimated error bars of this study. The 
estimated errors would have to be much larger for the models to be in agreement. For 
reference, the result for the CMB model fitting two resonances with all open channels is 
shown in Table lU Note that the errors quoted for the CMB results include systematic 
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1.8 1 1.2 

Energy (GeV) 

FIG. 3: Real parts (left panels) and imaginary parts (right panels) of vrA^ ttN (upper panels) 
and ttN r]N (lower panels) scattering amplitudes in the Sn partial wave. Dotted lines give the 
resonant, dashed lines the non-resonant, and solid lines the total scattering amplitude. Partial- 
wave amplitudes are shown with error bars. The resonant model used here is the two-resonance, 
two-channel CMB model described in the text, with non-resonant amplitudes described by distant 
poles. 



errors from model uncertainties, and so are larger than the errors arising only from fitting 
uncertainties in the results of the present work. Many resonance parameters found in the 
two-channel CMB model are within the estimated errors of the full model. The exception 
is the width of the 5*11(1650), where the full CMB model value is larger than that of any of 
the other eight models. The full CMB model fit apparently has a complicated interference 
between the second resonance and the channels that are excluded in the truncated versions. 
We also list the recommended physical parameters for these states from the Particle Data 
Group . Almost all results are within the conservative estimated ranges they publish. It is 
interesting that although the results for the 5*11 (1535) full width have significant variations 
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1.8 1 1.2 1.4 

Energy (GeV) 

FIG. 4: Caption as in Fig.|31 except the resonant model used here is the X-matrix model described 
in the text, with non-resonant amplitudes described by distant poles. 



using the different models studied here, the values are all in the lower part of the PDG 
estimated range. Similarly, the CMB and fT-matrix results for i^^^r are all in the lower part 
of the PDG estimated range for that quantity. 

Breit-Wigner models with relativistic and non-relativistic resonance shapes are fitted 
separately and give almost identical results. Despite having the best fits, they give results 
that vary the most from the PDG values. Compared to the other results, the Breit-Wigner 
models have significantly lower mass, and smaller widths and vrA^ branching fractions for 
S'ii(1650). This is more true for the fits using the polynomial background than for those 
using the distant poles non-resonant amplitudes. Since the resonant peak in Im(T) smoothly 
blends in with the non-resonant background, the S'ii(1650) width is very sensitive to how 
this background is treated. Therefore, it is not surprising that the biggest discrepancies 
arise in this resonance property. Interference with the overlapping lower energy state also 
has an important influence. 
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FIG. 5: Caption as in Fig. 13 except the resonant model used here is the Breit-Wigner model with 
non-relativistic widths described in the text, with non-resonant amplitudes described by distant 
poles. 



Comparisons between the CMB and i^-matrix models are the most interesting. The full 
fit amplitudes in Figs. El and |3] are very similar. They even miss the sharp structure in 
the ttN elastic amplitude at ~ 1.6 GeV in the same way. At first glance the resonant 
and non-resonant amplitudes are also very similar, but the small difference in the shape 
of the imaginary elastic amplitude is the primary source of the differences in the extracted 
full width in the two models. The S'ii(1535) full width from the CMB model with distant 
poles background is larger than for the i^'-matrix model with distant poles, and vice versa 
for the S'ii(1650) full width. For the error bars we derive, the difference in the full width 
is more significant than the difference in Bt^m. This shows the interplay of the interfering 
resonances. On the other hand, the same two models with polynomial background (where 
the CMB model no longer satisfies analyticity) have very similar values for the full width. 
Interpretation of this result will be a key component of the discussion. 
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FIG. 6: Caption as in Fig. 13 except the resonant model used here is the Breit-Wigner model with 
relativistic widths described in the text, with non-resonant amplitudes described by distant poles. 



IV. DISCUSSION 



A variety of empirical models were used to fit partial wave amplitudes (the input 'data') 
for a carefully chosen problem, the Su partial wave and its two most important channels, 
ttN and t]N. This partial wave is interesting because of the physics interest in the Sn reso- 
nances and the large uncertainty in their properties as reported by the PDG ^|. By using 
identical input amplitudes and fitting strategies, we have made the first objective compari- 
son of A^* resonance extraction models. Four different resonance models (CMB, i^-matrix, 
and Breit-Wigner with non- relativistic and relativistic widths), and two different empirical 
models for the non-resonant amplitude (distant poles and polynomial) are employed. These 
models are in regular use for the extraction of hadronic properties. It should be emphasized 
that the four resonance models have almost identical amplitudes for isolated resonances. 
The primary differences among the models come from the way the dynamics of resonance 
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Energy (GeV) 

FIG. 7: Caption as in Fig.|31 except the resonant model used here is the two-resonance, two channel 
CMB model described in the text, with non-resonant amplitudes described by polynomials. 

interference, multichannel effects, and non-resonant amplitudes are treated. Although res- 
onance models of widely varying quality are employed, the non-resonant models are both 
rather empirical. However, this partial wave in these channels has non-resonant amplitudes 
which are comparatively small. Since the non-resonant amplitudes have a smooth energy 
dependence, their influence on the extracted resonance parameters should be small. The 
main purpose of this work is to study the model dependence of the extraction of Sn prop- 
erties (mass, width, and vrA^ branching fraction) in a case where overlapping resonances, 
multichannel effects, and analyticity constraints are all expected to be important. 

The models used here can be put in order according to the theoretical constraints em- 
ployed. The CMB model with distant poles background satisfies multichannel unitarity and 
constraints from analyticity, and handles resonance-resonance quantum mechanical interfer- 
ence well. In fact, the CMB model includes the most complete resonance propagation effects 
of any of the existing models. The i^-matrix model used here does not satisfy analyticity 
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Energy (GeV) 

FIG. 8: Caption as in Fig.|ll except the resonant model used here is the X-matrix model described 
in the text, with non-resonant amplitudes described by polynomials. 



constraints, and leaves out re-scattering dynamics present in the CMB model. The Breit- 
Wigner models used here satisfy essentially no theoretical constraints. However, unitary 
Breit-Wigner models 3] and i^-matrix models that satisfy analyticity constraints have 
been developed. 

The large range of properties for the two lowest energy Su resonances in the Review of 
Particle Properties [l(l| is also found here. This is evidence that much of the uncertainty in 
the PDG estimates of 5*11 properties comes from model dependence, since the same input 
amplitudes are used in every fit. No evidence is found for a third Su state in the energy 
range studied. Some of the Breit-Wigner models have the best fits to the data, but this 
is due to the flexibility of these models rather than an ability to describe the underlying 
dynamics. Arbitrary adjustments must be employed in order to obtain good flts to the data. 
The empirical phases between the resonances provide a simple way to adjust the resonance- 
resonance interference at the cost of obscuring the physics output. As a result, the physical 
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FIG. 9: Caption as in Fig. |S1 except the resonant model used here is the Breit-Wigner model 
with non-relativistic widths described in the text, with non-resonant amplitudes described by 
polynomials. 

properties of the S'ii(1650) determined with the Breit-Wigner models are very different than 
with the other models. 

The models with the strongest theoretical constraints, the CMB and i^-matrix models 
with distant poles background, provide better agreement with each other and with the CMB 
fits to a much larger set of reactions One major result of the present work is the 

differences between CMB and i^-matrix models found in a situation where their differences 
should be the largest. The extent of the disagreement cannot be simply stated. For the 
error bars given for the input amplitudes and those determined in our fits, the difference 
between the S'ii(1535) mass and full width for the CMB and i^-matrix models with distant 
poles is significant. This is an important measure of the difference in these two models. 

Other considerations could contribute to these differences. The estimated errors quoted 
in the partial- wave amplitude fits have a direct effect on the values quoted here. It is possible 
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FIG. 10: Caption as in Fig.|31 except the resonant model used here is the Breit-Wigner model with 
relativistic widths described in the text, with non-resonant amplitudes described by polynomials. 



that the errors in the partial wave amplitudes are understated due to a lack of understanding 
of the model dependence. Although a more sophisticated non-resonant amplitude could be 
required, the smoothness of these amplitudes in all models argues against this. Each model 
has problems that are likely due to the truncated channel space employed. The most obvious 
problem is with the width of 5*11(1650). About 20% of the overall strength that was supposed 
to go to channels other than vrA^ and riN has to be included somewhere in the smaller channel 
space. 

In the context of the small but important test case chosen for this study, the results of 
the CMB and /^-matrix model fits are found to have small, but potentially important dif- 
ferences. Since the CMB model is better constrained theoretically, the resonance properties 
extracted using this model should be preferred when there is reason to doubt other models. 
A possible conclusion could be that this proves that the constraints provided by analyticity, 
a superior treatment of re-scattering, and an improved treatment of resonance interference 
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are important in this case. Prom the present hmited comparison, such a conclusion is hkely 
premature. This issue requires further study, particularly from the theoretical side. 

We should note that many N* states will not be obscured by strong threshold effects 
and will not have strong interference with other states of the same angular momentum and 
parity. The simplified dynamics of the /^-matrix model then give it the practical advantage 
of a simpler and more stable path to a good fit to the partial-wave amplitudes. Proper use of 
this model in the analysis of N* data is unhkely to give results that are badly incorrect. Por 
resonances with a larger number of open channels, the interaction of the non-resonant and 
resonant amplitudes can be much more complicated and simplified fitting can give erroneous 
results. 

The primary result of this paper is that even in a small multi-channel problem, dynamics 
are important. Since Breit-Wigner models have very few theoretical constraints, ad-hoc 
parameters are required to fit real data such as those of the two-channel problem studied 
here. We therefore suggest that consideration of these issues be part of any attempt to 
determine global recommendations for baryon resonance properties. 
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